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How the X-ray spectrum of a Narrow-Line Seyfert 1 galaxy may be 
reflection dominated 
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ABSTRACT 

A model for the inner regions of accretion flows is presented where, due to disc instabili- 
ties, cold and dense material is clumped into deep sheets or rings. Surrounding these density 
enhancements is hot, tenuous gas where coronal dissipation processes occur. We expect this 
situation to be most relevant when the accretion rate is close to Eddington and the disc is 
radiation-pressure dominated, and so may apply to Narrow-Line Seyfert 1 (NLS1) galaxies. 
In this scenario, the hard X-ray source is obscured for most observers, and so the detected 
X-ray emission would be dominated by reflection off the walls of the sheets. A simple Comp- 
tonization calculation shows that the large photon-indices characteristic of NLSls would be 
a natural outcome of two reprocessors closely surrounding the hard X-ray source. We test 
this model by fitting the XMM-Newton spectrum of the NLS1 1H 0707-495 between 0.5 and 
1 1 keV with reflection dominated ionized disc models. A very good fit is found with three 
different reflectors each subject to the same Y = 2.35 power-law. An iron overabundance is 
still required to fit the sharp drop in the spectrum at around 7 keV. We note that even a small 
corrugation of the accretion disc may result in Y > 2 and a strong reflection component in 
the observed spectrum. Therefore, this model may also explain the strength and the variabil- 
ity characteristics of the MCG-6-30-15 Fe Ka line. The idea needs to be tested with further 
broadband XMM-Newton observations of NLS Is. 
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1 INTRODUCTION 

The basic model for the active nucleus of a typical Seyfert 1 galaxy 
consists of a black hole surrounded by an ultraviolet and soft X- 
ray emitting thin accretion disc above which a patchy active corona 
emits hard X-rays. The model is supported by strong X-ray reflec- 
tion signatures due to Compton backscattering and fluorescence of 
the harder coronal X-ray s by the dense disk ( Pounds et al. 1990[ 
Nandra & Pounds 1994; Petrucci et al. 200C). In some objects an 



thought to explain the unusual properties of Narrow-Li ne Seyfert 1 
(NLS1) galaxies (e.g., founds, Done & Osborne 1995 ). 



Recently, Boiler et al. (2002 ) reported on a "sharp spectral fea- 



extreme red wing is found in the iron line of the reflection spec- 



trum (Ianaka et al. 1995; Nandra et al. 1999; Fabian et al. 2000), 



showing that the disk can extend very close to the black hole where 
large relativistic effects occur. In general the strength of the reflec- 
tion spectrum indicates that the disk surface is approximately flat, 
although disk instabilities could well make it ribbed or clumpy, par- 
ticularly if the disk^sj^idiation^rjre^ 

& Eardley 1 974; |Guilbert & Rees 19881; |Krolik 1998;; T urner, Stone 
& Sano 2001). Such a situation is most likely to occur when the ac- 
cretion rate is high and close to the Eddington limit, which is often 
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ture" at about 7.1 keV in the XMM-Newton spectrum of the NLS1 
1H 0707-495 (z = 0.0411). This deep drop (over a factor of two 
in flux) in the spectrum occurs at almost the exact energy of the 
neutral iron edge, so absorption models (in particular, a partial cov- 
ering model) were favoured to explain the feature. However, these 
models suffered from a very soft intrinsic power-law (T ~ 3.5) and 
an unreasonably large value of the iron abundance (~ 35 x solar). 
Moreover, it was difficult for the partial covering model to account 
for the absence of a narrow Fe Ka line and to explain the rapid 
variability of the source over the entire waveband. The sharp drop 
in the spectrum could also be due to the blue wing of a relativistic 
iron emission line, but, as the authors noted, to explain the depth 
of the drop at ~7 keV requires "invoking a very extreme Fe abun- 
dance and/or reflection fraction." Here, we consider the possibility 
of a reflection dominated X-ray spectrum in detail. 

The basic idea we wish to explore is that within a system 
which is accreting at close to the Eddington rate, disc instabili- 
ties funnel denser material into many deep rings, between which 
the hard X-rays are emitted by some coronal process (Figure [l]). To 
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Figure 1. The proposed geometry. The inner region of the radiation- 
pressure dominated accretion disc is broken up into a series of deep sheets 
or rings, with a typical size, R, of the order of the disc thickness, and 
separated by a distance d. Hard X-ray sources are found in between the 
sheets, but are obscured from the view of observers along most line-of- 
sights. Therefore, the detected X-ray emission is dominated solely by the 
reflected emission from the walls of the sheets. In Section^, we have con- 
sidered circular sheets of radius R instead of squared ones to simplify the 
computations. 



an off-axis observer, the X-ray emission will consist mostly of a 
reflection component. The rapid variability seen in the X-ray emis- 
sion from 1H 0707^-95 could then in part be due to changes in 
the geometric obscuration of one sheet on another, ex plaining why 
little spectral variability is seen ( jBoller et al. 2002| ). This reflec- 
tion model requires a less extreme input spectrum and hidden flux 
than the partial-covering model. While this is clearly an idealized 
picture, recent numerical simulations of the magneto-rotational in- 
stability in radiation-pressure dominated discs have found that the 
turbulence can cause density contrasts up to a factor of 200 (T urner, 
Stone & Sano 2001). Therefore, it does seem plausible that a 
radiation-pressure dominated disc will exhibit strong density inho- 
mogeneities. 

In the next section, we calculate what the hard X-ray power- 
law may look like in our idealized picture and show that a steep 
spectrum naturally arises from this situation. In Section ^| we fit 
reflection dominated spectra to the broadband XMM-Newton data 
of 1H 0707^-95. Finally, we discuss the implications of this model 
in Section ^. 



2 CALCULATION OF THE X-RAY POWER-LAW 

Here we give a simple analytic estimate of the slope of the X-ray 
continuum by assuming that in the inner radiation-pressure domi- 
nated part of the accretion disc the dense and cold sheet-like part 
of the flow i s embedded in a tenuous, hot plasma (cf., Guilbert 
& Rees 1 988; |Celotti, Fabian & Rees 1992|; Kuncic, Blackman & 
Rees 1996 ^Kuncic7 Celotti & Rees 1997| )! The accretion energy 
is dissipated, probably via magnetic reconnection, in the hot phase 
only, with the cold sheets reprocessing (reflecting and thermaliz- 
ing) the incident hard flux. 

For ease of computation, the picture in Fig. [I] was adapted so 
that the cold sheets have a flat disc-like geometry with radius R. 



We consider two such discs facing each other and separated by a 
distance d. The hot phase is sandwiched by the two cold discs and 
therefore has a cylindrical slab geometry. From a point in the mid- 
plane of the cylinder, at a distance r from its centre, the covering 
fraction of the two cold sheets is 
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The average total covering fraction is then given by 
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With the hard luminosity of the hot phase being Lh, the re- 
processed luminosity is given by 



L Tcp — CLh(1 — I 



(3) 



where a is the albedo of the cold sheets. We make use of the re- 



sults of Vlalzac, Beloborodov & Poutanen (2001 ) who used Monte 
Carlo simulations of Comptonizing coronae above dense cold ma- 
terial and showed that the albedo depends on the spectral index of 
the illuminating radiation, being smaller for steep spectra. The ther- 
mal reprocessed soft photon flux from the cold sheets is assumed to 
be a black-body spectrum with temperature r rcD , an d is calculated 
by the method outlined by Merloni & Fabian (2001). The hard lu- 
minosity is assumed to be due only to inverse Compton scattered 
photons and can be ap proximated as a sum of a cut- off power-law 
and Wien component (|Wardzinski & Zdziarski 2000). 



The emerging spectral index of the power-law is given by 
lnP sc 



1 = a = 



ln(l + 49 + 166 2 



(4) 



where & = kT e /m e c 2 is the electron temperature of the hot phase, 
and P sc is the scattering probability averaged over the source vol- 
ume and depends only on the corona l optical depth r. In a slab 
geometry P sc can be approximated as (Zdziarski 1994) 



Psc = 1 + 



exp(- 



t J It 2 w ' 



(5) 



where E\ is the exponential integral. 

Once the geometry is fixed (by fixing the aspect ratio R/d), 
the temperature of the hot phase, and consequently the slope of 
the Comptonized continuum, can be calculated self-consistently by 
solving Lh = Lc = f^Z an Lc(x)dx. We have calculated the 
spectral index T for different values of the ratio R/d and in Figure □ 
we plot the spectral index for three different values of the optical 
depth of the hot phase. Clearly, as the covering fraction increases, 
the spectrum softens significantly. Also, the spectrum is steeper for 
a denser hot medium (larger r). Finally, if there is accretion energy 
dissipated in the cold phase, t hen the spectrum will be softer for a 
given R/d (cf., jvlalzac 200 1| ). 

Figure |3| shows an equivalent plot where the abscissa is now 
the half opening angle between the two reprocessors 9/2. These 
plots show that spectral indices typical of NLSls (i.e., T > 2.1- 
2.2) are a natural consequence of this scenario if the half opening 
angle between the two reprocessors is less than about 20 degrees. In 
such a case, it is likely that a observer may only detect the reflected 
emission arising from the walls of the reprocessors. 
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Figure 2. The spectral index T as a function of the aspect ratio R/d, for 
r = 0.1 (solid line), r = 1 (dashed line) and r = 10 (dot-dashed line). 
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Figure 4. The upper panel shows the folded three-reflector model and 
EPIC-pn data of 1H 0707-495 between 0.5 and 1 1 keV. The lower panel 
displays the residuals to the fit in units of standard deviations. 
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Figure 3. The spectral index T as a function of the half opening angle of 
the couple of disc-like cold sheets for r = 0.1 (solid line), r = 1 (dashed 
line) and r = 10 (dot-dashed line). 



3 APPLICATION TO THE NLS 1 1H 0707^*95 

In this section, we attempt to fit the broadband X-ray spectrum of 
the NLS1 1H 0707-495 with models of reflection dominated spec- 
tra, which might be expected if the accretion disc has broken into a 
number of dense segments (Fig. |l|). We consider the EPIC-pn spec- 
trum of 1H 0707^195 between 0.5 and 1 1 keV and fit wit h the con- 
stant density ionized disc mo dels of |Ross & Fabian (1993 ) (see also 
Ross, Fabian & Young 19991. The parameters of the model are the 



ionization parameter £ = 4?tJ 7 x/wh, where Fx is the illuminating 
flux between 10 eV and 100 keV, the photon index of the power-law 
spectrum that strikes the slab, and a normalization constant. Rela- 
tivistic bl urring was a pplied to the model during fitting using the 
kernel of Laor (1991). The emissivity index was fixed at —3, but 



the inner and outer radii, along with the inclination angle, were left 
free. Absorption due to the Galactic column of 5.79 x 10 20 cm -2 
was included in all fits. 

If one tries to fit the data with a single reflection spectrum, 
then the best fit (x 2 /d.o.f.=506/309) is obtained with a model with 
a three times solar Fe abundance. Recall that no power-law compo- 
nent has been added to these models, giving them, in effect, an in- 
finite reflection fraction. We find that T = 2.44 and log £ = 2.040, 
which results in a strong neutral Fe Ka line. Extreme Kerr blurring 
(rtn = 1-24 r B ) was needed to fit the large red wing of the line be- 
tween 2.5 and 6 keV. However, as indicated by the large x 2 , there 
were significant residuals to the fit, particularly below 1 keV where 



the model p redicts strong emission lines from Fe and O (R oss & 
Fabian 1993). To dilute the effect of these lines, a power-law was 
added to the model with the value of V fixed to be the same as 
the reflected component. Reflection still dominated the spectrum 
with a reflection fraction > 10. The power-law improved the best 
fit to x 2 /d.o.f.=479/308, but required an Fe abundance greater than 
5 times solar. In this fit, T softened slightly to 2.52, but there was 
little change in the ionization parameter, the inner and outer radii, 
or the inclination angle (~ 23 degrees). There are still residuals 
below 1 keV, but, interestingly, there is now a clear line-like resid- 
ual just above 6 keV. Adding a Gaussian to the fit drops the x by 
38 with the addition of three degrees of freedom - significant at 
> 99.99 per cent, according to the F-test. Here, the best-fit model 
has a seven times overabundance of Fe, and T = 2.56. The line had 
an energy of 6.62 keV and a width a — 0.235 keV, indicating that 
it probably also arises from somewhere in or on the disk. 

To investigate this further, the power-law and Gaussian com- 
ponents of the model were replace by a second reflector. The second 
reflection component was drawn from the same group of models as 
the first, and was also subject to relativistic blurring. We fixed the 
photon-index and inclination angle to be the same for both com- 
ponents. The addition of this second reflector greatly improved the 
fit (x 2 /d.o.f.=342/305), as the combined emission removed most of 
the residuals from the soft X-ray band, and now only a five times 
overabundance of Fe was required. The ionization parameter of the 
second component is log£ = 3.397 resulting in a strong He-like 
Fe Ka line at 6.7 keV and an inner radius of 5.4 r g . The other re- 
flector remained relatively neutral with log£ = 2.09 and strongly 
blurred (n n = 1.24 r g ). The photon-index flattened slightly to 
r = 2.41. While this fit is almost acceptable there is still a line-like 
residual above 6 keV. Adding a Gaussian component to account 
for this feature significantly improves the fit (x 2 /d.o.f.=292/302). 
This new line has an energy of 6.45 keV and a width of 0.358 keV, 
demonstrating that the line arises from weakly ionized Fe in the 
accretion disc. 

The two-reflector plus Gaussian emission line model gives a 
very good fit to the EPIC-pn data between 0.5 and 1 1 keV, how- 
ever it is interesting to consider replacing the line component with 
a third reflection spectrum. This model also results in an excel- 
lent fit to the data (x 2 /d.o.f. =287/303; Fi gure ft), but requires a 7 
times solar abundance of Fe. The three components all have dif- 
ferent ionization parameters: log£ — 1.986 (with ri n = 2.08 r g 
and r out = 7.25 r g ), 3.847 (with r m = 4.18 r g ), and 1.722 (with 
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Figure 5. The three-reflector model that was successfully fit to 1H 0707— 
495. The solid line shows the total model, the dashed line plots the highly 
blurred but neutral inner reflector, the dot-dashed line shows the ionized 
reflector with r ln = 4.81 r g , and the dotted line plots the outer neutral 
reflector with r ln = 19.8 r g . These models do not contain any contribution 
from the incident power-law with T = 2.35. 



Tin = 19.8 r g ). The last two reflectors had their outer radii frozen 
at 50 r g . The photon index and the inclination angle for all three 
components were fixed to have the same values. The best fit values 
were V = 2.35 and 18.2 degrees, respectively. The model is shown 
in Figure H. 



4 DISCUSSION 

We have shown that the X-ray spectrum of 1H 0707-495 is well 
fitted by a reflection-dominated spectrum. The incident power-law 
is mildly steep with a photon index of 2.35. Furthermore, the model 
can completely account for the strong soft excess in this object by 
a combination of blurred emission lines and thermal emission from 
the irradiated blobs. The geometry envisaged for the source consists 
of a set of deep dense sheets or rings, orbiting the central black hole, 
between which hard X-rays are emitted. The fact that at least two 
reflectors are needed to provide a good fit to the data is not incon- 
sistent with the model. There should be numerous rings or sheets 
in the centre of the disk with (as in the standard coronal picture) 
multiple hard X-ray sources active at any one time. The system is 
observed at sufficient inclination that the hard X-ray sources are 
not directly visible and only X-rays reflected from the sides of the 
dense matter are seen. 

In fact, the model we have presented is a simplified picture 
of an intrinsically complex situation. In this respect, the interpreta- 
tion of the fit parameters should be consider with some caution. In 
a more realistic situation, the reflection features will be broadened 
due to both relativistic effects and disc turbulence, so that parame- 
ters like ri n and r out should not be taken at face value, but may be 
interpreted as the signature of the presence of different reproces- 
sors at different locations. Also, a range of ionization parameters is 
to be expected from the complex configuration we envisage, as a 
reflector would likely have an ionization gradient on its surface. 

A consistent model for both the photon index of the input 
power-law and the reflection component is obtained if the open- 
ing angle of the gaps between the sheets or rings is less than 40 
degrees and the whole disk inclination is about 20 degrees (Fig. ^ 
and Sect. ^). The metallicity is still required to be high (about 5-7 



times the Solar value) but less than for other models. It is possible 
that the metallicity can be lowered if multiple reflection is included 
(this will have an approximately similar effect to squaring the re- 
flection spectrum). 

The spectral-independent variability of 1H 0707^-95 (B oiler 
et al. 2002) would be mostly due to changes in the brightness of the 
hard X-ray source causing the observed reflection components to 
vary. Additional variability may also arise from changes in the ge- 
ometry of the system, with sheets or rings obscuring each other as 
they orbit the black hole. The geometry may also account for larger 
scale variations. When observed by XMM-Newton, 1H 0707-495 
was ab out ten times fainter than in previous observations (B oiler et 
al. 2002). This could be due to slight changes in the geometry and 
not necessarily a large change in intrinsic luminosity. If the hard 
X-ray source was previously more directly visible (say the opening 
angle of the sheets or rings was larger, or the hard X-ray emitting 
corona extended nearer to the disk surface) then the observed flux 
would have been larger. Indeed, the s pectrum obser ved by ASCA 



when the flux was higher (r ~ 2.27, Leighly 1999), is consistent 
with being the unobscured continuum emission peeking through 
the cold sheets. 

It is possible that the model we have developed here for 
1H 0707-495 has a wider relevance to Seyfert 1 galaxies. The re- 
flection fraction can exceed unity if the disk is in the form of deep 
rings or sheets, and leads to an intrinsic hard power-law steeper than 
r = 2. MCG-6-30-15 has sometimes been observed with a spec- 
tral index of 2.1 or steeper (|Vaughan & Edelson 200 1^ Shih, Fabian 
& Iwasawa 2001) and i ts strong iron line suffers f ar less variabil- 



ity than its continuum (Vaughan & Edelson 2001). This could be 



due to corrugations in its inner disk which may cause the reflec- 
tion components to be stronger and more long-lasting than if the 
accretion disc was perfectly flat. 
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